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ABSTRACT 

Detection of gamma-ray bursts (GRBs) from redshifts z J> 7 would open a new window into the 
earliest epoch of cosmic star formation. We construct separate star formation histories at high redshifts 
for normal (Pop I and II) stars, and for predominantly massive (Pop III) stars. Based on these separate 
histories, we predict the GRB redshift distribution to be observed by the Swift mission. Regardless 
of whether Pop III progenitors are able to trigger GRBs, we find that a fraction ~ 10% of all bursts 
detected by Swift will originate at z ^ 5. This baseline contribution is due to Pop l/ll star formation 
which must have extended out to high redshifts in rare massive galaxies that were enriched by heavy 
elements earlier than the typical galaxies. In addition, we consider the possible contribution of Pop III 
progenitors to the observable GRB rate. Pop III stars are viable progenitors for long-duration GRBs 
which are triggered by the collapsar mechanism, as long as they can lose their outer envelope through 
mass transfer to a companion star in a close binary. We find that the likelihood of Pop III binaries to 
satisfy the conditions required by the collapsar mechanism could be enhanced significantly relative to 
Pop l/ll binaries. If Pop III binaries arc common, Swift will be the first observatory to probe Pop III 
star formation at redshifts z > 7. 

Subject headings: binaries: general — cosmology: theory — gamma rays: bursts — stars: formation 



1. INTRODUCTION 

The first stars in the universe, so-called Population III 
(hereafter Pop III), formed out of metal- free gas at the 
end of the cosmic dark ages or redshifts z J> 10 (for re- 
views, see, e.g., Barkana & Loeb 2001; Bromm & Larson 
2004; Ciardi & Ferrara 2005; Glover 2005). These stars 
are predicted to have been predominantly very massive 
with M* £ 100Af Q (Bromm, Coppi, & Larson 1999, 2002; 
Abel, Bryan, & Norman 2002; Nakamura & Umemura 
2001), and to have left a mark on the thermal and chem- 
ical evolution of the intergalactic medium (IGM). First, 
their predicted high surface temperatures (e.g., Bond, Ar- 
nett, & Carr 1984) implies that they may have been effi- 
cient sources of ionizing photons (e.g., Tumlinson & Shull 
2000; Bromm, Kudritzki, & Loeb 2001b; Schaerer 2002). 
A contribution from Pop III stars to the reionization of 
the IGM may be required to account for the large optical 
depth to Thomson scattering inferred by the Wilkinson 
Microwave Anisotropy Probe (WMAP; Kogut et al. 2003) 
during its first year of operation (e.g., Cen 2003; Wyithe & 
Loeb 2003a, 2003b). Second, since the stellar evolutionary 
timescale for massive Pop III stars is short ~ 10 6 yr, the re- 
sulting initial enrichment of the IGM with heavy elements 
could have occurred rather promptly (e.g., Mori, Ferrara, 
& Madau 2002; Bromm, Yoshida, & Hernquist 2003; Wada 
& Venkatesan 2003; Yoshida, Bromm, & Hernquist 2004). 

Gamma-Ray Bursts (GRBs) offer unique prospects for 
probing the cosmic star formation (Totani 1997; Wijers 
et al. 1998; Blain & Natarajan 2000; Porciani & Madau 
2001; Bromm & Loeb 2002; Hernquist & Springel 2003; 
Mesinger, Perna, & Haiman 2005; Natarajan et al. 2005) 
as well as the IGM (Loeb 2003; Barkana & Loeb 2004; 



Gou et al. 2004; Inoue, Omukai, & Ciardi 2005; Ioka 
& Meszaros 2005) at redshifts z > 7, beyond the cur- 
rent horizon of galaxy and quasar surveys. GRBs are the 
brightest electromagnetic explosions in the universe (for 
recent reviews, see, van Paradijs, Kouveliotou, & Wijers 
2000; Meszaros 2002; Piran 2004) and their emission is de- 
tectable out to z > 10. The detectability of their gamma- 
ray emission (Lamb & Reichart 2000) allows simultaneous 
monitoring of a major portion of the sky, while the de- 
tectability of their afterglow (Ciardi & Loeb 2000) allows 
to determine their high redshift from the appearance of the 
intergalactic Lya absorption trough in the infrared (Loeb 
2003; Barkana & Loeb 2004). The popular collapsar model 
for the central engine of long-duration GRBs (MacFadyen, 
Woosley, & Heger 2001 and references therein), involves 
the collapse of a massive star to a black hole (BH). This 
model naturally explains the observed association of long- 
duration GRBs with star-forming regions (e.g., Fruchter 
et al. 1999; Djorgovski et al. 2001; Bloom, Kulkarni, & 
Djorgovski 2002a) and the Type Ib/c supernova signature 
on the spectra of rapidly decaying afterglow lightcurves 
(e.g., Bloom et al. 2002b; Hjorth et al. 2003; Matheson et 
al. 2003; Stanek et al. 2003). 

Because of their high characteristic masses, Pop III stars 
could potentially lead to high redshift GRBs. The re- 
cently launched Swift satellite 3 (Gehrels et al. 2004) is 
ideally suited to utilize this novel window into the high- 
rcdshift universe. In the following sections we address the 
underlying question: which fraction of high- redshift bursts 
could originate from Pop III progenitors ? The actual frac- 
tion and distribution of high- z GRBs to be measured by 
Swift, might reflect the absence or presence of the potential 



Pop III contribution. 
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The organization of the paper is as follows. In § 2 we 
describe our cosmic star formation model, in particular de- 
termining the Pop III mode at the highest redshifts. The 
resulting GRB redshift distribution is calculated in § 3, 
together with a discussion of plausible GRB progenitors. 
Finally, we address the implications of our results in § 4. 

2. STAR FORMATION MODES AT HIGH REDSHIFTS 

We note that GRBs are expected to exist at redshifts 
z J> 7 even in the absence of any true Pop III contribu- 
tion. This is due to the rapid enrichment of the IGM with 
heavy elements dispersed by the first supernovae (SNe) be- 
ginning at z ;> 20 (e.g., Loeb & Haiman 1997; Madau, Fer- 
rara, & Rees 2001; Bromm et al. 2003; Furlanetto & Loeb 
2003; Scannapieco, Schneider, & Ferrara 2003). Once a 
given region of the universe has been enriched beyond a 
critical metallicity, Z alt <~ 10~ 3 ' 5 Z Q , the mode of star for- 
mation is predicted to shift from high-mass Pop III stars 
to the lower-mass Pop I and II cases (e.g., Omukai 2000; 
Bromm et al. 2001a; Schneider et al. 2002; Bromm & 
Loeb 2003b; Mackey, Bromm, & Hernquist 2003; Schnei- 
der et al. 2003). The mass fraction of super-critical gas 
is rapidly growing toward lower redshift, and GRBs could 
be formed in the conventional way from metal-enriched 
Pop I and II progenitors at high-z. In the following dis- 
cussion, we first construct the total cosmic star formation 
rate (SFR), and subsequently decompose the total SFR 
into separate Pop I/II and Pop III components. 

2.1. Star Formation History 

Our model for the total cosmic SFR closely follows that 
of Bromm & Loeb (2002), and we here only briefly de- 
scribe the key assumptions. The abundance and merger 
history of the cold dark matter (CDM) halos is described 
by the extended Press-Schechter formalism (Lacey & Cole 
1993). We assume that the IGM has a two-phase structure, 
consisting of neutral and ionized hydrogen phases. The 
reionization of the IGM was likely an extended process, 
occurring over 6 <; z <; 20 (e.g., Cen 2003; Wyithe & Loeb 
2003a; Sokasian et al. 2004; Furlanetto & Loeb 2005). To 
bracket the possibilities, we consider two reionization red- 
shifts, z ro i on i=a 7 and 17, where z ro ; on corresponds to an 
ionization filling fraction by volume of ~ 50%. In each 
case, reionization is spread out over a range in redshifts, 
Az/(l + z)~ 1. 

Within each phase of the IGM, stars are able to form in 
two different ways. The first mechanism pertains to pri- 
mordial, metal-free, gas. Such gas undergoes star forma- 
tion provided that it accretes onto a dark matter halo with 
a sufficiently deep gravitational potential well or equiva- 
lently a mass above a minimum value. For the neutral 
medium, this minimum mass is set by the requirement that 
the gas will be able to cool. Radiative cooling by molecular 
hydrogen (H 2 ) allows star formation in halos with a virial 
temperature T v ; r J> 500 K, while atomic cooling domi- 
nates for halos with T v j r ;> 10 3 9 K. Since H 2 can be easily 
photo-dissociated by photons below the Lyman-limit, its 
significance in the cosmic star formation history is unclear 
(e.g. Bromm & Larson 2004 and references therein), and 
so we only show results without H 2 cooling in this paper. 
We note, however, that in the limiting case of negligible 
H 2 photodissociation feedback, the cosmic star formation 



rate at z ;> 15 could be larger by one order of magnitude 
than the purely atomic cooling case discussed here. 
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Fig. 1. — Cosmic comoving star formation rate (SFR) in units of 
Mq yr _1 Mpc~ 3 , as a function of redshift. We assume that cool- 
ing in primordial gas is due to atomic hydrogen only, and the star 
formation efficiency is r) t = 10%. (a) Late reionization (2 rc ion ~ 7). 
Solid line: Total comoving SFR. Dotted lines: Contribution to the 
total SFR from Pop I/II and Pop III for the case of weak chem- 
ical feedback. Dashed lines: Contribution to the total SFR from 
Pop I/II and Pop III for the case of strong chemical feedback, (b) 
Early reionization (z rc ion ~ 17). We adopt the same convention 
for the lines as in panel (a). In all cases, Pop III star formation 
is restricted to high redshifts, but extends over a significant range, 
Az ~ 10 - 15. 



For the ionized medium, on the other hand, the min- 
imum threshold mass is given by the Jeans mass, since 
the infall of gas and the subsequent formation of stars re- 
quires that the gravitational force of the dark matter halo 
be greater than the opposing pressure force on the gas. 
After reionization, the IGM is photo-heated to tempera- 
tures ;> 10 4 K, leading to a dramatic increase in the Jeans 
mass. We model the suppression of gas infall according 
to results from spherically-symmetric collapse simulations 
(see Bromm & Loeb 2002 for details). In calculating the 
late reionization case (z ro i on w 7), we employ the predic- 
tion by Thoul & Weinberg (1996) that gas infall is com- 
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pletely suppressed in halos with circular velocities v c <^ 
35 km s _1 . For the early reionization case (z ro ion ~ 17), 
however, we use the recent work by Dijkstra et al. (2004), 
showing that the infall suppression due to photo-ionization 
heating could be much less severe in the high-redshift uni- 
verse. Specifically, we assume that in this latter case infall 
is completely suppressed only in halos with circular veloc- 
ities v c <; 10 km s _1 . 

Within our model, the second mechanism to form stars 
occurs in gas that has experienced a previous burst of 
star formation, and is therefore already somewhat enriched 
with heavy elements. Such gas, residing in a halo of mass 
Mi, can undergo induced star formation triggered by a 
merger with a sufficiently massive companion halo of mass 
M2 > 0.5Mi. We finally assume that stars form with an 
efficiency of 77* ~ 10%, independent of redshift and regard- 
less of whether the gas is primordial or pre-enriched. This 
efficiency yields roughly the correct fraction of fie found 
in stars in the present-day universe. Figure 1 shows the 
resulting total star formation histories. It is evident that 
there are two distinct epochs of cosmic star formation, one 
at z <~ 3, and a second one at z ~ 8 for late reionization, 
whereas there is only a single, extended peak at z <~ 5 for 
early reionization. 

2.2. Population III Star Formation 

To determine the fraction of the total SFR contributed 
by Pop III stars, we have to identify those halos that cross 
the atomic cooling threshold for the first time. In addi- 
tion, we require that the collapse takes place in a region 
of the IGM which is not yet enriched with heavy elements 
from previous episodes of star formation. Here, we adopt 
the formalism developed in Furlanetto & Loeb (2005) who 
derived the redshift-dependent probability that a newly 
collapsed halo forms out of pristine gas (see their Fig. 2). 
This probability crucially depends on the efficiency with 
which the newly created metals are dispersed into the IGM 
via SN driven winds. To bracket the range of possibilities, 
we consider the cases of weak and strong chemical feed- 
back, corresponding to winds experiencing large and small 
radiative losses, respectively 4 . 

As can be seen in Figure 1, Pop III star formation is 
limited to the highest redshifts, but in each case extends 
over a substantial range in redshift: Az ~ 10 for strong 
chemical feedback, and Az ~ 15 for weak feedback. The 
Pop III histories are rather similar for both early and late 
reionization. The suppression of gas infall for the early 
reionization case (with v c <J 10 km s _1 ), would have a 
much more pronounced effect on halos that cool via H 2 , 
because of their shallower potential wells. 

3. GRB REDSHIFT DISTRIBUTION 

Next we will predict the GRB redshift distribution for 
flux-limited surveys, distinguishing between the contribu- 
tions from Pop I/II and Pop III star formation. In partic- 
ular, we will focus on the existing Swift satellite, which is 
capable of making the most detailed determination of the 
GRB redshift distribution to date. 

3.1. Population I/II Contribution 

4 The weak and strong feedback cases correspond to K^f 3 = 1/3 and 



Assuming that the formation of GRBs follows closely 
the cosmic star formation history with no cosmologically 
significant time delay (e.g, Conselice et al. 2005), we write 
for the number of all GRB events per comoving volume 
per time, regardless of whether they are observed or not: 
V'grbM = ^grb x il>*(z), where ip*(z) is the cosmic SFR, 
as calculated in §2. The efficiency factor, ?7gr.b, links the 
formation of stars to that of GRBs, and is in principle a 
function of redshift as well as the properties of the under- 
lying stellar population. The stellar initial mass function 
(IMF) is predicted to differ fundamentally for Pop I/II 
and Pop III (e.g., Bromm & Larson 2004 and references 
therein). The GRB efficiency factor will depend on the 
fraction of stars able to form BHs, and consequently on 
the IMF (see §3.2). Here, we assume that ?7grb is constant 
with redshift for Pop I/II star formation, whereas Pop III 
stars may be characterized by a different efficiency. 

The number of bursts detected by any given instrument 
depends on the instrument-specific flux sensitivity thresh- 
old and on the poorly-determined isotropic-luminosity 
function (LF) of GRBs (see, e.g., Schmidt 2001; Sethi & 
Bhargavi 2001; Norris 2002). In order to ascertain what 
Swift is expected to find, we modify the true GRB event 
rate as follows: 



VgrbW=%rb^W / p(L)dL . (1) 

J Llim(z) 



Here, p(L) is the GRB LF with L being the intrinsic, 
isotropic-equivalent photon luminosity (in units of pho- 
tons s _1 ). If /ii m denotes the sensitivity threshold of a 
given instrument (in photons s _1 cm -2 ), then the mini- 
mum luminosity is: 



L lim (z) = 4nd 2 L / Um (l + z) a - 2 , (2) 



where <1l is the luminosity distance to a source at redshift 
z, and a the intrinsic high-energy spectral index (Band et 
al. 1993). For definiteness, we assume a = 2, which gives 
a reasonable fit to the observed burst spectra (see Band 
et al. 1993 for a detailed discussion). We here use the 
same lognormal LF and the same parameters as described 
in Bromm & Loeb (2002). 

In Bromm & Loeb (2002), we predicted that ~ 25% of 
all bursts observed by Swift would originate at z > 5. This 
estimate was based on a flux threshold of /ii m = 0.04 pho- 
tons s _1 cm~ 2 . Based on the first few months of actual ob- 
servations by Swift, the sensitivity limit has recently been 
revised upward to fu m — 0.2 photons s _1 cm~ 2 , compa- 
rable to the older BATSE experiment (e.g., Berger et al. 
2005). Using this revised flux limit, we predict that <~ 10% 
of all Swift GRBs will originate at z > 5. 

1 in equation (4) of Furlanetto & Loeb (2005). 
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Fig. 2. — Predicted GRB rate observed by Swift. Shown is the ob- 
served number of bursts per year, dN9^ B /dhi(l + z), as a function 
of redshift. All rates are calculated with a constant GRB efficiency, 
VGRB — 2 X 10 — 9 bursts Mg , using the cosmic SFRs from Fig. 1. 
(a) Late reionization (z r cion ~ 7). Dotted lines: Contribution to 
the observed GRB rate from Pop I /II and Pop III for the case of 
weak chemical feedback. Dashed lines: Contribution to the GRB 
rate from Pop I /II and Pop III for the case of strong chemical feed- 
back, (b) Early reionization (z r cion ~ 17). The lines have the same 
meaning as in panel (a). Filled circle: GRB rate from Pop III stars 
if these were responsible for rcionizing the universe at z ~ 17 (see 
text). The GRB rates from Pop III progenitors are very uncertain; 
they could be zero, or, in the other extreme, display an enhancement 
by up to one order of magnitude above the baseline rates shown here 
(see text for details). 

Over a particular time interval Ai b s in the observer 
frame, the observed number of GRBs originating between 
redshifts z and z + dz is 



dN GRB _ jobs / \ A *ob: 

Vgrb\ z ) 



dV 



dz 



(1 + z) dz 



(3) 



where dV/dz is the comoving volume element per unit 
redshift (see Bromm & Loeb 2002). As a final step, we 
normalize the GRB formation efficiency per unit mass in 
Pop I/E stars to t] GRB ~ 2 x 10~ 9 GRBs M" 1 . This 
choice results in a predicted number of <~ 90 GRBs per 



year detectable by Swift. In Figure 2, we show the Swift 
GRB rate, associated with Pop I/II star formation. For 
both early and late reionization, the observed distribution 
is expected to peak around z <~ 2. This distribution is 
broadly consistent with the first GRB redshifts, still lim- 
ited in number, measured during the first months of the 
Swift mission (Berger et al. 2005). We now turn to the 
possible contribution to the high redshift GRB rate from 
Pop III stars. 

3.2. Population III Contribution 

We begin by assuming that Pop III star formation gives 
rise to GRBs with the same (constant) efficiency as is em- 
pirically derived for Pop I/II stars. As is evident from Fig- 
ure 2, only for the case of weak chemical feedback is Swift 
expected to detect a few bursts deriving from Pop III pro- 
genitors over the ~ 5 yr lifetime of the mission. Whether 
reionization occurred early or late, on the other hand, has 
only a small effect on the predicted rates. It is quite pos- 
sible, on the other hand, that Swift will not detect any 
Pop III GRBs at all. Regardless of the uncertain contri- 
bution from Pop III stars, however, the prediction that 
~ 10% of all Swift bursts originate at z > 5 is rather ro- 
bust. This fraction is due to Pop I/II progenitors that are 
known to produce GRBs, and those should exist at z > 5. 

Adopting the same 7?grb for Pop III as for Pop I/II, 
however, could be significantly in error. To examine the 
fundamental difference between the stellar populations, we 
need to go beyond the phenomenological approach pursued 
so far and discuss the properties of plausible GRB progen- 
itors in greater physical detail. 

3.2.1. Collapsar Engine 

Existing evidence implies that long-duration bursts are 
related to the death of a massive star, leading to the for- 
mation of a BH (see review by Piran 2004). The pop- 
ular collapsar model assumes that an accretion torus is 
temporarily formed around the black hole, and that the 
gravitational energy released during the accretion is able 
to power a strong explosion (e.g., Woosley 1993; Lee & 
Ramirez-Ruiz 2005). For the collapsar to result in a GRB, 
additional requirements have to be met beyond the for- 
mation of a BH. We will discuss these next, and will then 
explore whether a subset of Pop III stars could successfully 
launch a GRB under the collapsar scenario. 

To successfully produce a GRB with a collapsar, three 
basic requirements have to be fulfilled (see, e.g., Zhang & 
Fryer 2004; Petrovic et al. 2005): 

(i) The progenitor star has to be sufficiently massive to 
result in the formation of a central BH. Collapse to a BH 
could occur either directly for initial masses of the pro- 
genitor ;> 4OM , or in a delayed fashion by fallback of 
the ejecta following a failed SN explosion for progenitor 
masses 25 £ M* <, 4OM (e.g., Heger et al. 2003). The 
number of BH forming stars resulting from a given total 
stellar mass, here denoted by ?7bh, will depend on the stel- 
lar IMF which in turn is predicted to differ between the 
Pop I/II and Pop III cases. 

For simplicity, we assume that the IMF in both cases 
consists of a power-law with the standard Salpeter value, 
dN/dm oc to -2 35 , but with different values for the lower 
and upper mass limits, M\ ovl and M up respectively. For 
Pop I/II stars, we take these to be: Mi ow = 0.1M Q and 
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M up = 1OOM . The Pop III IMF, on the other hand, 
is still very uncertain (see, e.g., Bromm & Larson 2004). 
The upper mass limit can be conservatively estimated to 
be M up ~ 500M o (Bromm & Loeb 2004), whereas for the 
lower limit, we consider two possibilities: M\ ow ~ 30M Q 
(e.g., Tan & McKee 2004) and - 100M Q (e.g., Abel et al. 
2002; Bromm et al. 2002). 
In general, 



M BH 



rn 



-2.35 



dm 



^ p m~ 1 - 35 dm 

"low 



(4) 



where M B h — 25M Q . For the Pop I/II case, this re- 
sults in ?7bh — 1/(700M Q ). The Pop III lower mass limit 
exceeds the threshold for BH formation in either case, 
Mi ow > Mbh- Not every Pop III star, however, will leave a 
BH behind. In the narrow mass range of ~ 140 — 260M Q , 
Pop III stars are predicted to undergo a pair-instability su- 
pernova (PISN) explosion (e.g., Fryer, Woosley, & Heger 
2001; Heger et al. 2003). A PISN will lead to the com- 
plete disruption of the star, such that no compact rem- 
nant will be left behind. For the Pop III case, the expres- 
sion results in tjbh — 1/(80M Q ) for M\ ow = 30M Q , and 
77bh ^ 1/(300M©) for M low = 100M Q . Thus, the BH for- 
mation efficiency is larger for Pop III compared to Pop I/II 
by up to one order of magnitude, depending on the lower 
mass limit. 

(ii) The progenitor star has to be able to lose its hydro- 
gen envelope in order for the relativistic outflow to pene- 
trate through and exit the star (e.g., Zhang et al. 2004). 
This requirement derives from the observed burst dura- 
tions, t <J 100 s, providing an estimate for the lifetime of 
the central GRB engine. The jet can therefore only travel 
a distance of r ~ ct ~ 5Oi?0 before being slowed down to 
non-relativistic speeds. Massive stars with hydrogen en- 
velopes grow to a large size during their later evolutionary 
phases. For example, red supergiants can reach radii of 
up to - 10 3 i? o (e.g., Kippenhahn & Weigert 1990). The 
effectiveness of mass loss crucially depends on metallicity 
(e.g., Kudritzki 2002), and on whether the star is isolated 
or part of a binary system. Below, we will discuss both 
effects further. 

(Hi) The progenitor star has to contain a central core with 
sufficient angular momentum to allow an accretion disk to 
form around the growing BH. Important aspects of stel- 
lar structure and evolution can be understood by dividing 
the star into a compact core, and an extended outer en- 
velope (e.g., Kippenhahn & Weigert 1990). Depending on 
the evolutionary stage, a radiative core is surrounded by a 
convective envelope, or vice versa. The pre-collapse stellar 
core has a mass M Cl radius R c , angular velocity w c , and is 
characterized by a specific angular momentum j c ~ R\u) c . 
Assuming that the collapse to a BH conserves specific an- 
gular momentum, the condition for an accretion torus to 
form around a growing BH in the center with mass Mbh is 
centrifugal support for material orbiting at the last-stable 
radius. This condition can be expressed as j c ;> j m i n , with 

j min = V6GM B h/c ~ 3 x 10 16 cm 2 s _1 (e.g., Podsiadlowski 
et al. 2004). 

Recent results obtained with sophisticated stellar evo- 
lution codes that include the effects of magnetic torques 
(e.g., Spruit 2002) have demonstrated the difficulty to 



identify progenitor systems for collapsar-driven GRBs 
that fulfill both conditions (ii) and (Hi). The basic prob- 
lem is that the removal of the extended H-envelope is ac- 
companied by the loss of angular momentum in the core 
(e.g., Petrovic et al. 2005). We will explore this problem 
next, first for single-star progenitors, and then for binaries. 

3.2.2. Single-star Progenitor 

In massive Pop I/II stars, radiation driven winds can 
lead to vigorous mass loss, where the main source of opac- 
ity is provided by metal lines (e.g., Kudritzki & Puis 
2000). Empirically, the existence of Wolf-Rayet (WR) 
stars proves that Pop I/II stars can indeed experience 
catastrophic mass loss, leading to the removal of the en- 
tire hydrogen-, and, in extreme cases, even of the helium- 
envelope. The violent mass loss, however, is accompanied 
by the effective removal of angular momentum from the 
remaining pre-collapse core, rendering the creation of a 
collapsar impossible. 

For massive Pop III stars, radiatively driven winds are 
predicted to be unimportant (e.g., Kudritzki 2002; Krticka 
& Kubat 2005). Alternatively, mass loss could occur as 
a result of stellar pulsations (through the e mechanism). 
Simplified, linear calculations, however, indicate that this 
mechanism is not important below ~ 500Mq (Baraffe, 
Heger, & Woosley 2001). There still remains an unex- 
plored possibility that Pop III stars could experience sig- 
nificant mass loss driven by radiation pressure on He + ions, 
where the opacity is provided by bound-free transitions. 
Alternatively, processed material from proceeding episodes 
of nuclear burning could be transported to the stellar sur- 
face by convection, rotationally-induced mixing, and dif- 
fusion, thus enriching the atmosphere to Z ;> 10~ 4 Z Q at 
which point line-driven mass loss is predicted to set in 
(Kudritzki 2002; Marigo, Chiosi, & Kudritzki 2003). Re- 
cently, it has been suggested that WR type winds, in con- 
nection with rapid rotation and the approach to the Ed- 
dington limit, could possibly lead to significant mass loss 
even for very low- metallicity stars (Vink & de Koter 2005). 
The most likely expectation, however, is that isolated mas- 
sive Pop III stars are not able to shed much mass prior to 
their final collapse. 

Thus it appears likely that the majority of massive, sin- 
gle star progenitors, both for Pop I/II and III, cannot give 
rise to a collapsar-driven GRB, although for different phys- 
ical reasons. We note, however, that the progenitor for 
the collapsar engine is still very uncertain. In addition to 
the binary scenario (see below), massive, rapidly rotating 
single stars have recently been considered (e.g., Yoon & 
Langer 2005; Woosley & Heger 2005). Our main results, 
based on the cosmic Pop III SFR and the IMF-dependent 
BH fraction, however, holds in this case as well. We next 
turn our attention to binary-star progenitors. 

3.2.3. Binary-star Progenitor 

Close binary systems provide a promising avenue to si- 
multaneously meet the requirements of strong mass loss 
combined with the retention of sufficient angular momen- 
tum in the collapsing core (e.g., Lee, Brown, & Wijers 
2002; Izzard, Ramirez-Ruiz, & Tout 2004). For Pop I/II, 
a binary pathway to a collapsar-driven GRB has already 
been suggested (e.g., Fryer et al. 1993). The basic idea is 
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that a close binary system, experiencing Roche-lobe over- 
flow (RLOF) when the primary evolves off the main se- 
quence, will go through a common-envelope (CE) phase, 
during which the hydrogen envelope of the primary can 
be removed without seriously draining away the spin of 
the remaining helium core (for a recent review, see Taam 
& Sandquist 2000). During the CE phase, the binary will 
spiral closer together, and the corresponding loss of orbital 
energy will heat the envelope with a given efficiency, often 
assumed to be very high: «ce — 1-0 (Taam & Sandquist 
2000). The frictional energy release during the in-spiral 
phase has been shown to be sufficient to unbind the hy- 
drogen envelope. 

For a progressively tightening binary, the spin of each 
component is tidally coupled to the orbital motion: lu c ~ 
w or b. Since the helium core is spun up again because of the 
spin-orbit coupling, it is able to retain sufficient angular 
momentum to fulfill the j c ;> j m j n requirement. 

3.2.4. GRB Formation Efficiency 
In general, we can express the GRB formation efficiency 
for Pop III stars as: 7y G RB ^ ^BHfybin^cioscfybcaming, where 
?7bin is the binary fraction, r) c \ ose the fraction of sufficiently 
close binaries to undergo RLOF, ^beaming — 1/50-1/500 
the beaming factor, where we conservatively assume that 
Pop III bursts are collimated by the BH central engine to 
the same angle as Pop I /II progenitors (the inferred colli- 
mation angles by Frail et al. 2001; Panaitescu & Kumar 
2001, are indeed comparable to those of radio jets from 
the much more massive BHs in galactic nuclei). We have 
already discussed ?7bh, and how Pop III star formation is 
characterized by an enhancement of up to one order of 
magnitude because of the higher fraction of BH-forming 
progenitors. 

It is currently not known whether Pop III stars can form 
binaries, and if so, what the corresponding binary frac- 
tion will be (e.g., Saigo, Matsumoto, & Umemura 2004). 
In present-day star formation, the incidence of binaries is 
high, with <~ 50% of all stars occurring in binaries or small 
multiple systems (e.g., Duquennoy & Mayor 1991). Cur- 
rent three-dimensional simulations of the formation of the 
first stars still lack the resolution to resolve the possible 
fragmentation of a collapsing cloud into tight binaries or 
multiple stellar cores on scales <J 100 AU. Although we 
cannot yet conclusively address the formation of close bi- 
naries, there is evidence from numerical simulations that 
binary or multiple clump formation is rather common on 
larger scales (;> 0.1 pc). In simulations where the collaps- 
ing gas had acquired a high degree of angular momentum, 
and where the collapse led to a disk-like configuration, 
pre-stellar clumps commonly occurred in binary or mul- 
tiple systems (e.g., Bromm et al. 1999, 2002; Bromm & 
Loeb 2003a, 2004). Thus motivated, our best guess is that 
f?bin <J 0.5. More work, in particular involving improved 
numerical simulations, is required to constrain this crucial 
quantity further. 

Provided that Pop III star formation does include a frac- 
tion of binaries, we can use the collapsar requirements to 
obtain an estimate for the maximum binary separation, 
&i,max, prior to the CE inspiral phase, as follows. Assum- 
ing for simplicity that a Pop III binary has equal mass 
components, the Roche radius is r^ ~ 0.5<Zj. A CE phase 
will only occur when the star overflows its Roche lobe 



during the red giant phase: Rrg > r L o>i- We esti- 
mate the Pop III radius during the giant phase, which is 
smaller than that of a Pop I star of equal mass, to be 
-Rrg ~ 300i?MS, where Rms is the main-sequence radius. 
Massive Pop III stars obey a simple mass-radius relation 
(Bromm et al. 2001b): M cx R^s- The maximum sep- 
aration for RLOF and therefore a CE phase to occur is 
thus 

/ M \ 1/2 

The minimum possible binary separation, ai^ min , on the 
other hand, will determine the extent of the inspiral pro- 
cess and therefore of the accompanying frictional heating 
of the hydrogen envelope as well as the tidal spinning up 
of the helium core. We approximately assume that the 
minimum separation is given by twice the radius of a mas- 
sive Pop III star during the main-sequence phase (which 
is only weakly dependent on mass): aj im i n ~ 10i? Q . If we 
further assume that the initial separations are distributed 
with equal probability per logarithmic separation interval, 
as they are for Pop I/II binaries (e.g., Abt 1983; Hea- 
cox 1998; Larson 2003), dN/d\na cx const., and that the 
largest possible binary separation is given by the Jeans 
length for the typical conditions in a primordial gas cloud, 
Xj ~ 1 pc (e.g., Bromm et al. 2002, Bromm & Loeb 2004), 
we estimate the fraction of all Pop III binaries that are suf- 
ficiently close to experience a CE phase to be ?7 c i oso ~ 30%. 

In summary, 77grb for Pop III stars is very uncertain, 
and could be zero in case that no Pop III binaries existed. 
On the other hand, one may argue that the binary prop- 
erties for Pop III had been similar to Pop I/II, in cases 
where Pop III star formation takes place in more massive 
host systems that could give rise to a stellar cluster, or at 
least multiple stars. Such a clustered environment is often 
suggested to explain the formation and the properties of 
binaries in present-day galaxies (e.g., Larson 2003). As- 
suming that the Pop III binary properties are similar to 
Pop I/II, we find a significant enhancement in ?7bh due to 
the difference in the IMF between the populations. We can 
then place an upper limit on the GRB rate from Pop III 
stars by multiplying the baseline rates in Figure 2 by a 
factor of <~ 10. This would result in Pop III GRB rates 
as large as ~ 10 bursts detected by Swift per year for 
the case of weak chemical feedback. Such very high rates 
can already be excluded, since Swift has only identified 
two GRBs'from z £ 5 as of yet, GRB 050814 at z ~ 5.3 
(Jakobsson et al. 2005) and GRB 050904 atz~ 6.3 (An- 
tonelli et al. 2005; Haislip et al. 2005; Kawai et al. 2005). 
For strong chemical feedback, on the other hand, we pre- 
dict rates of less than one burst detected per year, even 
if the BH efficiency were increased by one order of magni- 
tude, and such a Pop III contribution cannot be excluded 
with the current constraints from Swift. 

3.3. Constraints from Reionization 

If massive Pop III stars led to an early reionization of 
the universe at z <~ 17, as may be required by the WMAP 
data (Kogut et al. 2003), we can obtain an estimate for the 
corresponding Pop III SFR at z ro j on and for the possible 
accompanying GRB rate, as follows. 

Pop III stars with masses ^ 100M Q produce <~ 10 62 H- 
ionizing photons per solar mass over their ~ 2x 10 6 yr 
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lifetime (e.g., Bromm et al. 2001b). We can then esti- 
mate that ~ 3 x 10 5 M Q in Pop III stars are required to 
produce ~ 5 ionizing photons for every hydrogen atom 
in a comoving Mpc 3 . This number is sufficient to com- 
pensate for recombinations at the mean cosmic density. 
Assuming further that the burst of Pop III star forma- 
tion is spread over a fraction e of the Hubble time at 
z rcion ~ 17, Aig F ~ 4x 10 7 (e/0.2) yr, the comoving 
Pop III SFR able to reionize the universe at that redshift 
is: SFR roion - 7.5 x 10- 3 /(e/0.2) M yr^ 1 Mpc" 3 . The 
extremely rapid growth of the collapsed fraction of baryons 
with redshift implies a value of e <C 1; however, the mini- 
mum value of e is ~ 0.1 because even within a single dark 
matter halo, star formation cannot be synchronized to bet- 
ter than the dynamical time which amounts to e ~ 0.1 for 
a virial density contrast of ~ 200. 

In Figure 2, we show the GRB rate which would corre- 
spond to SFR rc i on for e = 0.2 when the constant Pop I/II 
GRB efficiency factor is used, resulting in ~ 0.1 GRBs 
per year. Under this conservative assumption, Swift is 
not expected to detect, within its expected 5 year mission 
lifetime, any bursts connected to the Pop III stars that 
were responsible for an early reionization of the universe. 
Again, the prospects for detection would be significantly 
improved if the Pop III GRB efficiency is boosted due to 
the increased fraction of BH-forming progenitors. 

4. SUMMARY AND CONCLUSIONS 

Figure 2 leads to the robust expectation that ~ 10% of 
all Swift bursts should originate at z J> 5. This prediction 
is based on the likely existence of Pop I/II stars in galaxies 
that were already metal-enriched at these high redshifts. 
Additional GRBs could be triggered by Pop III stars, with 
a highly uncertain efficiency. Assuming that long-duration 
GRBs are produced by the collapsar mechanism, a Pop III 
star with a close binary companion provides a plausible 
GRB progenitor. We have estimated the Pop III GRB 
efficiency, reflecting the probability of forming sufficiently 
close and massive binary systems, to lie between zero (if 



tight Pop III binaries do not exist) and ~ 10 times the em- 
pirically inferred value for Pop I/II (due to the increased 
fraction of BH-forming progenitors among Pop III stars). 

Recently, Gorosabel et al. (2004) and Natarajan et al. 
(2005) predicted the expected redshift distribution of long- 
duration GRBs, assuming they trace the cosmic star for- 
mation history, with various phenomenological prescrip- 
tions for the dependence on metallicity. In difference from 
these studies, we isolate the zero-metallicity (Pop III) stars 
and treat them as potential GRB progenitors based on a 
physical model. 

It is of great importance to constrain the Pop III star 
formation mode, and in particular to determine down to 
which redshift it continues to be prominent. The extent 
of the Pop III star formation will affect models of the ini- 
tial stages of reionization (e.g., Wyithe & Loeb 2003a, b; 
Ciardi, Ferrara, & White 2003; Sokasian et al. 2004; 
Yoshida et al. 2004; Alvarez, Bromm, & Shapiro 2006) 
and metal enrichment (e.g., Mackey et al. 2003; Furlan- 
etto & Loeb 2003, 2005; Scannapieco et al. 2003; Schaye 
et al. 2003; Simcoe, Sargent, & Rauch 2004), and will 
determine whether planned surveys will be able to effec- 
tively probe Pop III stars (e.g., Scannapieco et al. 2005). 
The constraints on Pop III star formation will also de- 
termine whether the first stars could have contributed a 
significant fraction to the cosmic near-IR background (e.g., 
Santos, Bromm, & Kamionkowski 2002; Salvaterra & Fer- 
rara 2003; Dwek, Arendt, & Krennrich 2005; Kashlinsky 
2005; Madau & Silk 2005; Fernandez & Komatsu 2006). If 
Pop III binaries were common, Swift might be the first in- 
strument to detect Pop III stars from galaxies at redshifts 
z >. 7. 
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